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We explore in detail the possibility of generation of continuous-variable (CV) entangled states 
of light fields with well-localized phases. We show that such quantum states, called entangled 
self-phase locked states, can be generated in nondegenerate optical parametric oscillator (NOPO) 
based on a type-II phase-matched down-conversion combined with polarization mixer in a cavity. 
A quantum theory of this device, recently realized in the experiment, is developed for both sub- 
threshold and above-threshold operational regimes. We show that the system providing relative 
phase coherence between two generated modes also exhibits different types of quantum correlations 
between photon numbers and phases of these modes. We quantify the entanglement as two-mode 
squeezing and show that CV entanglement is realized for both unitary, non-dissipative dynamics 
and for dissipative NOPO in the entire range of pump field intensities. 
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INTRODUCTION 

It is now believed that entanglement of quantum com- 
posite systems with continuous degree of freedom is the 
basis of most applications in the field of Quantum In- 
formation 0]. Interest in continuous variable (CV) en- 
tanglement is being extensively excited by successful ex- 
periments on quantum teleportation based on two-mode 
squeezed states as well as the experiments dealing with 
entanglement in atomic ensembles 0. Since then, re- 
markable theoretical and experimental efforts have been 
devoted to generating and quantifying CV entangled 
states. 

In this paper we propose a novel type of CV en- 
tangled states of light-field, which contain well-localized 
phases. They arc different from the well-known entan- 
gled Einstein-Podolsky- Rosen (EPR) states generated in 
a nondegenerate optical amplifier jj, 0, which exhibit 
large phase fluctuations. We believe that such entangled 
states can be generated in a self-phase locked nondegener- 
ate optical parametric oscillator (NOPO), based on the 
type-II phase-matched down-conversion and additional 
phase localizing mechanisms stipulated by the intracav- 
ity waveplate. The motivations for this study are the 
following: 

For the first time the CV entangled states of light 
were studied in and demonstrated experimentally 
in for nondegenerate optical parametric amplifier 
(NOPA). Then a CV entanglement source was built from 
two single-mode squeezed vacuum states combined on 
a beamsplitter 0. It is well known that each of the 
orthogonally polarized and frequency degenerate fields 
generated by NOPO is a field of zero-mean values. The 
phase sum of generated modes is fixed by the phase of 
the pump laser, while their phase difference undergoes 
a phase diffusion process 0] stipulated by vacuum fluc- 
tuations. As a rule, the NOPO phase diffusion noise is 



substantially greater than the shot noise level, that limits 
the usage of NOPO in precision phase-sensitive measure- 
ments. Various methods based on phase locking mecha- 
nisms 00 00| have been proposed for reducing such 
phase diffusion. In the comparatively simple scheme real- 
ized in the experiment 0] , self-phase locking was achieved 
in NOPO by adding an intracavity quarter-wave plate to 
provide polarization mixing between polarized signal and 
idler fields. The evidence of self-phase locking was pro- 
vided there by the high level of phase coherence between 
the signal and idler fields. Following this experiment, 
the semiclassical theory of such NOPO was developed 
in |7J. Recently, the schemes of multiphoton parametric 
oscillators based on cascaded down-conversion processes 
in x' 2 ' media placed inside the same cavity and show- 
ing self-phase locking have been proposed 0. As was 
demonstrated in , the system based on combination of 
OPO and second harmonic generation also displays self- 
phase locking. The formation of self-phase locking and 
its connection with squeezing in the parametric four-wave 
mixing under two laser fields has been demonstrated in 
[Io| . An important characteristic of self-phase locked de- 
vices concerns the phase structure of generated subhar- 
monics. Indeed, the formation of the variety of distinct 
phase states under self-phase locked conditions has been 
obtained in the mentioned Refs. 0000 0] . It was 
recently noted that the schemes involving phase locking 
are potentially useful for precise interferometric measure- 
ments and optical frequency division because they com- 
bine fine tuning capability and stability of type - II phase 
matching with effective suppression of phase noise. That 
is why we believe it will be interesting to consider phase 
locked dynamics also from the perspective of Quantum 
Optics and, in particular, from the standpoint of produc- 
tion of CV entanglement. 

A further motivation for such task is connected with 
the problem of experimental generation of bright entan- 
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gled light. So far, to the best of our knowledge, there 
is no experimental demonstration of CV entanglement 
above the threshold of NOPO. The progress in experi- 
mental study of bright two-mode entangled state from 
cw nondegenerate optical parametric amplifier has been 
made in |llj. The theoretical investigation of CV entan- 
gled light in transition through the generation threshold 
of NOPO is given in [l2j • One of the principal experimen- 
tal difficulties in advance toward a high-intensity level is 
the impossibility to control the frequency degeneration 
of modes above the threshold. We hope that the usage 
of phase locked NOPO may open a new interesting pos- 
sibility to avoid this difficulty. 

In this paper we report what is believed to be the first 
investigation of self- phase locked CV entangled states. 
We develop the quantum theory of self-phase locked 
NOPO, with decoherence included, in application to 
the generation of such entangled states. This scheme 
is based on the combination of two processes, namely, 
type-II parametric down-conversion and linear polariza- 
tion mixer with cavity-induced feedback. The paramet- 
ric down-conversion is a standard technique used to pro- 
duce an entangled photon pairs as well as CV two-mode 
squeezed states @j- The beam splitter including polar- 
ization mixer are also considered as experimentally ac- 
cessible devices for production of entangled light-fields 
Besides these, there have been some studies of a 
beam splitter for various nonclassical input states, in- 
cluding two-mode squeezing states |l4j . It is obvious, 
and also follows from the results of the mentioned pa- 
pers 0, 0, , that the operational regimes of the com- 
bined system with cavity-induced feedback and dissipa- 
tion drastically differ from those for pure processes. We 
show below that analogous situation takes place in the 
investigation of quantum-statistical properties of a com- 
bined system such as the self-phase locked NOPO. 

The paper is arranged as follows. In Section II we 
formulate the model of combined NOPO based on the 
processes of two-photon splitting and polarization mix- 
ing, and present a semiclassical analysis of the system. 
Section III is devoted to the analysis of quantum fluctua- 
tions of both modes within the framework of linearization 
procedure around the stable steady-state. In Section IV 
we investigate the CV entangling resources of self-phase 
locked NOPO on the base of two- mode squeezing for both 
sub-threshold and above-threshold operational regimes. 
We summarize our results in Section V. 



MODEL OF SELF-PHASE LOCKED NOPO 

As an entangler we consider the combination of two 
processes in a triply resonant cavity, namely, the type - 
II parametric down-conversion in medium and po- 
larization mixing between subharmonics in lossless sym- 
metric quarter- wave plate. The Hamiltonian describing 



intracavity interactions is 

H = ihE (e i( -* L -"V4 - e-^-^aa) 
+ihk (e'* fc a^a^ — e'~ l ^ k a^a\a 2 ) 
+h X {e^a+a 2 +e^a ia +) , (1) 

where eij are the boson operators for the cavity modes 
LUi. The mode 03 at frequency to is driven by an external 
field with amplitude E and phase while aiand a 2 
describe subharmonics of two orthogonal polarizations 
at degenerate frequencies lu/2. The constant fce** fe de- 
termines the efficiency of the down-conversion process. 
Linear coupling constant denoted as x e ** x describes the 
energy exchange between the modes and besides, \ l& de- 
termined by the amount of polarization rotation due to 
the intracavity waveplate, <!> x determines the phase dif- 
ference between the transformed modes. We take into 
account the detunings of subharmonics Aj and the cav- 
ity damping rates 7$ and consider the case of high cavity 
losses for pump mode (73 3> 7). However, in our analy- 
sis we take into account the pump depletion effects. The 
reduced density operator p within the framework of the 
rotating wave approximation and in the interaction pic- 
ture is governed by the master equation 

dp 1 

-qT = tt [H inU p] + ^2 7i (2a t pa+ - a+a lP - pafa,) 

i=l,2 

k 2 

H (2aia 2 pa,ia 2 — afaia 2 a 2 p — paf a\a 2 a 2 ) , (2) 

73 

where 

Hi n t = hAia^ai + hA 2 a 2 a 2 
kE 

+ ih {ai a 2 ~ a i a 2) + hx{ a i a 2 + a\a 2 )■ (3) 

73 

Let us also note that this equation is rewritten through 
the transformed boson operators dj —> ciiexp(— 
with being $3 = $ L , $2 = \ (*L + $fc - $1 = 
\ ($l + &k + $ x )- That leads to cancellation of the 
phases on the intermediate stages of calculations. As a 
result, the Hamiltonian Hint depends only on real- valued 
coupling constants. The appearance of the last term in 
the master equation means that the influence of the adi- 
abatically eliminated fundamental mode is reduced to an 
additional loss mechanism for the subharmonic modes. 

In order to proceed further, we now consider the phase- 
space symmetry properties of the system. It is easy to 
check that the interaction Hamiltonian satisfies the com- 
mutation relation [Hi n t, U (ir)] — with the operator 
U (9) = exp \i0 (afai — 0.^0,2)] . Moreover, analogous 
symmetry [p (t) , U (tt)} = takes place for the density 
operator of the system, which obeys the master equation 
(0- Using this symmetry we establish the following se- 
lection rules for the normal-ordered moments of the mode 
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operators: 



{a{ aj a\ai) = 0, 



(4) 



if k + I + m + n ^ 2j, j = 0, ±1, ±2, .... Another pe- 
culiarity of NOPO with Hamiltonian (JIJ is displayed 
in the phase space of each of the subharmonic modes. 
The reduced density operator for each of the modes 
is constructed from the density operator p by tracing 
over the other mode Pi(2) = ^ r 2(i) (p)- Therefore, we 
find that [ Pl (t) , Ui {n)] = [p 2 (t) , U 2 {■*)} = 0, where 
Ui (9) = exp (i9a~lai) , (i — 1,2) is the rotation opera- 
tor. It is easy to check using these equations that the 
Wigner functions W\ and W 2 of the modes have a two- 
fold symmetry under the rotation of the phase-space by 
angle 7r around its origin, 



W t (r, 9) = W t (r, 9 + tt) , 



(5) 



where r, 9 are the polar coordinates of the complex phase 
space. Note, that such symmetry relationships (J5J 
radically differ from those taking place for the usual 
NOPO without the quarter-wave plate. Indeed, in this 
case, (case of x = in Eqs. 0-©) the Wigner functions 
Wi are rotationally symmetric and the symmetry prop- 
erties © read as (a^a 1 ^™ 1 ^) = 0, if k - I ^ m - n. 

We perform concrete calculations in the positive P- 
representation \n\ in the frame of stochastic equations 
for the complex c-number variables a{ and Pi correspond- 
ing to the operators ai and af 



da! 

~dt = " (71 



«Ai) ai 
+ (e— \a1a2) 



^ = _( 7l -iA!)A 
+ (e - \P1P2) a 2 



i\a 2 + Ri, 



- ixfo + Rt 



(G) 



(7) 



Here e = kE/jz, X — fc 2 /73- The equations for a 2 , (3 2 
are obtained from ©, Q by exchanging the subscripts 
(1)<^(2). Rx t 2 are Gaussian noise terms with zero means 
and the following nonzero correlators: 

(R 1 (t)R 2 (t')) = (e- \ ai a 2 ) 5 (t - t 1 ) , (8) 



{R+{t)R+{t')) = (e- 6 (t - . (9) 

In this approach the stochastic amplitude is given by 
0:3 = (E — kaiot2) /j3- So, Eqs.©-© involve the deple- 
tion effect of the pump mode, which leads to the appear- 
ance of the above-threshold operational regime. 

First, we shall study the steady-state solution of the 
stochastic equations in semiclassical treatment, ignoring 
the noise terms for the mean photon numbers rijo and 



phases ipjo of the modes (rij — ctjftj, <fj = ^ In (aj/(3j)). 
The mean photon numbers read as 



"10 



""20 



1 /A2 

A Ui 

1 (Ax 
A V A 2 



1/2 r 



1/2 r 



{etrf + 1 2 



+ 7 2 — 7 



where 



7i ( A 2 
2 Ui 



1/2 



72 / Ai 
2 \A 2 

,± 



1/2 



(10) 



(11) 



(12) 



As we can see from (tTUl) . ijTT)) . 

positive for e exceeding two critical points 



10 and n 20 are real and 



(sty 



7172 + AiA 2 + x 



TV4x 2 AiA 2 - ( 7l A 2 - 72A 1 )^ ) (13) 



and, besides, the solutions n^ and n^,, (i=l,2) corre- 
spond to two distinct values e~£ r and e~ r accordingly. The 
steady-state values of the phases corresponding to each 
of the critical points e~ r , e^ r are obtained as 



sin(^20 - fto) = sin (^20 - fw) = 



1 

2^ 



7i 



A2 
Ai 



1/2 



■72 



Ai 



1/2 



cos(ipf Q + ipf ) = ^e 2 - (e'e 



+ 7 2 



(14) 



(15) 



It is easy to check that these solutions exist for both 
modes only if the following relation holds 



4x 2 AiA 2 > (71 A 2 - 72 Ai) 



(16) 



Let us note that the steady-state solutions (|14H . I|15|) com- 
pletely determine the absolute phases of the orthogonally 
polarized modes, which are hence self-locked, unlike the 
ordinary NOPO. These results are in accordance with 
the ones obtained in [(J , but for another configuration 
of NOPO. In the scheme proposed in only the signal 
and idler modes are excited in the cavity, while the pump 
field is a travelling wave. Nevertheless, in the adiabatic 
regime considered here there is correspondence between 
both schemes. Indeed, it is not difficult to check that the 
results ljTU|) . ljTT|) transform to the corresponding results 
of the mentioned scheme Q by replacing the parameter 
e with the corresponding pump field amplitude. 

We now turn to the standard linear stability analysis of 
these solutions, assuming for simplicity, the perfect sym- 
metry between the modes, provided that the cavity decay 
rates and the detunings do not depend on the polariza- 
tion ( 71 = 72 = 7, Ai = A 2 = A). The stability of the 
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system is governed by the matrices F and F + , F_ describ- 
ing the dynamics of small deviations 5at and 5 Pi from 
the semiclassical steady-state solutions (see, Sec. III). We 
reach stability if the real parts of eigenvalues of these 
matrices are positive. This analysis displays an evident 
dependence on the sign of the detunings Ai = A2 = A. 
For the positive detuning, A > 0, only the steady-state 
solutions n± Q — and y>J — ¥>i~0' + V10 are stable, 
while for the case of negative detuning the stability holds 
for the solutions with the (+) superscript. As this anal- 
ysis shows, for either sign of detuning, the threshold is 
reached at e ^ £th, where 



e th = ^(x-IAlf+T 2 , 



(17) 



and the steady-state stable solution for mean photon 
numbers can be written in the general form as 



1 

no = n w = n 2 Q = - 



e 2 - ( X - |A| 



(18) 



The phases are found to be 



V10 = ¥20 = -^Arcsin - (\ + |A|) + irk, (19) 

for A > 0. For the opposite sign of the detuning, A < 0, 
as we noted, the mean photon numbers are given by the 
same Ea. (|18|l . while the phases read as 

<Pio = ^csini(x+|A|)+7r^fc + , (20) 

1P20 = ^Arcsia^ (%+ |A|) +n (k - ^ , 

(k = 0, 1, 2, ..). In the region e ^ eth the stability condi- 
tion is fulfilled only for the zero amplitude steady-state 
solution a% = a 2 = fli = [3 2 = 0. So, the set of above- 
threshold stable solutions for both modes have two-fold 
symmetry in the phase-space which was indeed expected 
from symmetry arguments (JSJ). 

Let us now consider the output behavior of the sys- 
tem for the special scheme of generation, when the cou- 
plings of in- and out-fields occur at only one of the ring- 
cavity mirrors. Taking into account that only the fun- 
damental mode is coherently driven by the external field 
with (a™) = -E/73, while the subharmonic modes are 
initially in the vacuum state, we obtain for the mean 
photon numbers (in units of photon number per unit 
time) nf = E 2 /2-/ 3 , n° ut = 2 7i n j0 (i = l,2) and 
hence n° ut = nif* = 2jn . Accordingly, parametric 
oscillation can occur above the threshold pump power 
P t h = huj 3 E 2 h /2j3, where the threshold value of the 

pump field is equal to Eth — -jf y (x — |A|) 2 + 7 2 . 

We are now in a position to study quantum effects in 
self-phase locked NOPO and will state the main results 
of the paper concerning CV entanglement. 



ANALYSIS OF QUANTUM FLUCTUATIONS 

The aim of the present section is to study the quantum- 
statistical properties of self-phase locked NOPO in linear 
treatment of quantum fluctuations. Quantum analysis 
of the system using P- representation is standard. A 
detailed description of the method can be found in . 
We assume that the quantum fluctuations are sufficiently 
small so that Eqs.©, can be linearized around the 
stable semiclassical steady state o>i (t) = a" + <5a, (t), 
Pi (t) = fif + 8 Pi (t). This is appropriate for analyz- 
ing the quantum-statistical effects, namely CV entan- 
glement, for all operational regimes with the exception 
of the vicinity of threshold, where the level of quantum 
noise increases substantially. It should also be empha- 
sized at this stage that in the above-threshold regime of 
self-phase locked NOPO the steady-state phases of each 
of the modes are well-defined in contrast to what hap- 
pens in the case of ordinary NOPO, where phase diffusion 
takes place. According to this effect, the difference be- 
tween the phases, as well as each of the phases, can not 
be defined in the above-threshold regime of generation 
of the ordinary NOPO. On the whole, the well founded 
linearization procedure can not be applied for this case. 
Nevertheless, the linearization procedure and analysis of 
quantum fluctuations for ordinary NOPO become pos- 
sible due to the additional assumptions about temporal 
behavior of the difference between the phases of the gen- 
erated modes [lfij . 

We begin with consideration of below-threshold oper- 
ational regime, for E < E t h, where the equations lin- 
earized around the zero-amplitude solution can be writ- 
ten in the following matrix form 

^Sof = -F^SoT + R» (ct,t) , (21) 

where /i = 1,2,3,4 and Sa^ = (Sa±, 5a 2 , Sa 3 , Sa^) — 
{6ai,6a2,6Pi,8fo), R>* = (R u R 2 ,Rf,Rt)- The 4x4 
matrix F^ is written in the block form 



F = 

with 2x2 matrices 

7 + iS, ix 



A, B 
B* A* 



A 



IX, 



P I °' 1 

y + i5 J ' V 1, 



The noise correlators are determined as 

(IP (a, t) R v (a, t')) = (a) S (t - t') 
with the following diffusion matrix 

B, 



D = 



0, B* 



(22) 



(23) 



(24) 



(25) 
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First, we calculate the temporal correlation functions 
of the fluctuation operators. The expectation values of 
interest can be written as the integral 



(SaP {t)5a v (f)) = 




F T being the transposition of the matrix F, The integra- 
tion over dr can be performed using the following useful 
formula for operators DF T — FD, obtained by straight- 
forward calculation. As a consequence, we arrive at the 
expression De F * = e Ft D. Finally we obtain 

{Saf (t) 8a u (t 1 )) = -\ (F^e^-^D) , (27) 

I 



and hence for t = t' 



(Sa»(t)Sa»(t)) = --(F- 1 D)^. (28) 



This formula, however, is not very convenient for prac- 
tical calculations. Therefore, we rewrite the correlation 
functions of the quantum fluctuations 5ai, 5 Pi in a more 
simple form through the two-dimensional column vectors 
6a = (Sai, 8a2) T , 6/3 = (5Pi,5p2) T - Upon performing 
the calculation, we finally arrive at 



Sa (Sa) 



Sa (spy 



S 4 - 4A V 

e 2 

2 (S* 4 - 4A V 



-2*A, S 2 
S 2 , -2 X A 

S 2 , -2 X A 
-2xA, S 2 



x(s 2 - 

A(S 2 



2A 2 

-2 X 2 



A(S 2 

x(s 2 - 



-2 X 2 
2A 2 



(29) 



where S 2 is introduced as 

S 2 = 1 2 + X 2 + A 2 



(30) 



Note, that S 2 > in the below-threshold regime. 

As an application of these results the mean photon 
number in the below-threshold regime can be calculated 
as follows, 



£ 2 5 2 



ni 



Tl2 = 



2 (S 4 - 4A 2 x 2 



(31) 



Next we focus on the mode locked regime, for E > Eth, 
considering the Eqs. @, m terms of the fluctuations 
Srii (t) = m (t) — riio and 6<fi (t) = ipi (t) — tpiQ of photon 
number and phase variables. In this regime the dynamics 
described by the linearized equations of motion actually 
decouples into two independent dynamics for two groups 
of combinations Sn± = brii ± Sni, 5ip± — Sif2 ± Sipi. In 
fact, one has 



d ( Sn + 
dt \Sip + 



d_/5n- 
di\5<p- 



= -F, 



= -F- 



Sn+ 
Sip+ 

Sri- 
Sif- 

the 



Rn-[ 

R v -\ 

Rn- 

Rip- 

diffusion 



(32) 



The drift, F±, and 
(RiQRjV)) = DlpS(t-t'), = n 



(R n (t)R m (t')) 
respectively 

F, = 



= D { n JS(t-t'),(n,m 

2Xn , 4n esin((p 2 o + <Pio) 
0, 2( 7 + An ) 



(33) 
matrices, 

ri-,ip-) are 



(34) 



/ 27, 4n xsinn/A 
I A/no, 



£>(+)_/ 4n 7, - 2esin(^20 + fio) 
-2esin((y520 + fio), -■y/m 



(35) 



(36) 



Due to the decoupling between (+) and (— ) 
combinations of the modes, we conclude that the 
following temporal correlation functions are equal 
to zero, (5ip + (t) Sep- (t')) = (5n+ (t) <5n_ (*')) = 

(Sn± (t)5(p T (t 1 )) — 0, and hence ((Sipi) 2 ) — ((S(p2) 2 



(8nx) y = ((5712) J- The other correlation functions 

can be calculated in the same way as described above for 
the below-threshold regime. The temporal correlation 
functions are derived as 



(Sn± (t) 

W± (*) 



and hence 



(Sn± (t'),5 V± (*')) 



■-F± e ± I \D±, 



= -±F7'D, 



(37) 



(38) 



Performing the concrete calculations for each of the 
cases A > and A < 0, we arrive at the following results 
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' ( Sn A (Zn K, n A - 1 ( 4«o [7(7 + An ) + (x- I A |) 2 ] , -An (x- | A |)s%n(A) \ , . 

.Up J ^+'^+7 " 4An (7 + An ) ^ -2An (x- | A |) S i 5 n(A), -A 7 J ' ^ 

A - 1 f 4Tt "X(x- I A |), 2 X7S i 5 n(A) \ 

UJ ^-'^-)/ - IlAU { Zxisign(A), J-( 7 2 - | A | ( x - | A |) ) ' ( 40 ) 



We see that the considered system displays different 
types of quantum correlations in terms of the stochas- 
tic variables, namely, between photon-number sum and 
phase sum in the modes, as well as between photon- 
number difference and phase difference in the modes. 
These results radically differ from the ones taking place 
for usual NOPO, where the correlation between photon- 
number difference and phase sum in the modes is realized. 
The results obtained indicate the possibilities to produce 
entanglement with respect to the new types of quantum 
correlations. 



CV ENTANGLEMENT IN THE PRESENCE OF 
PHASE LOCKING 

Let us now turn our attention to quantum statistical 
effects and the entanglement production for the case of 
perfect symmetry between the modes (71 =72 = 7, Ai — 
A2 = A). We note that unlike the two-mode squeezed 
vacuum state, the state generated in the above-threshold 
regime of NOPO is non-Gaussian, i.e. its Wigner func- 
tion is non-Gaussian ^t|- Recently, it has been demon- 
srated that some systems involving beam splitters 
also generate non-Gaussian states. The general consider- 
ation of this problem for self-phase locked NOPO seems 
to be very complicated. However, the mentioned results 
allow us to conclude that the state generated in self-phase 
locked NOPO is most probably non-Gaussian. So far, the 
inseparability problem for bipartite non-Gaussian state 
is far from being understood. On the theoretical side, 
the necessary and sufficient conditions for the separa- 
bility of bipartite CV systems have been fully developed 
only for Gaussian states, which are completely character- 
ized by their first and second moments. To characterize 
the CV entanglement we address to both the insepara- 
bility and strong EPR entanglement criteria U| which 
could be quantified by analyzing the variances of the 
relevant distance V- = V (X\ — X2) and the total mo- 
mentum V+ = V (Yi + Y 2 ) of the quadrature amplitudes 
of two modes Xk — [a£ exp (— iOk) + dfe exp (i0h)\ > 

Yk = 73 \- a k exp - a k ex P (iOk)] ,{k=l, 2), where 

V(X) = (V 2 ) — (X) 2 is a denotation for the variance and 
9k is the phase of local oscillator for the k-th mode. The 
two quadratures and Yfe are non commuting observ- 
ables. The inseparability criterion for the quantum state 



of two optical modes reads as [l9j 

V=±(V+ + V-)<1, (41) 

i.e. indicates that the sum of variances drops below the 
level of vacuum fluctuations. Since the states of the sys- 
tem considered are non-Gaussian, the criterion (|41|l is 
only sufficient for inseparability. The strong EPR entan- 
glement criterion is quantified by the product of variances 
as V + V_ < j. We remind that the sufficient condition for 
inseparability l|41|) in terms of the product of variances 
reads as V+V- < 1, i.e. is weaker than the strong EPR 
condition. 

In order to obtain the general expressions for variances, 
we first write them in terms of the boson operators corre- 
sponding to the Hamiltonian . We perform the trans- 
formations ai — ► ai exp which restore the previous 
phase structure of the intracavity interaction. Using also 
the symmetry relationships (@J) we find quite generally 
the variances at some arbitrary quadrature phase angles 
0i, 9 2 as 

V± = V ±i?cos(A6>), (42) 

where 

V = -(V+ + V-) = l + 2n 

-2 I (aia 2 ) I cos(£0 + $ arg ), (43) 

R = 2Re((a 2 ) e iT,e ) - 2 | (a^a 2 ) |, (44) 

Ae = e 2 -e 1 - $ x , = e x + e 2 + * t + $ fe , (45) 

and n — (a^ai^ = ^a^a 2 ) is the mean photon num- 
ber of the modes, $ a rg = arg (0102). So, in accordance 
with the formula (|42|l , the relative phase <& x between the 
transformed modes gives the effect of the rotation of the 
quadrature amplitudes angle 62 — 9i ■ 

Obviously, the variances V± and hence the level of CV 
entanglement depend on all parameters of the system 
including the phases. The minimal possible level of V 
is realized for an appropriate selection of the phases 8i, 
namely for 0\ + 9 2 = — arg {aia 2 ) — — <&k- Further, in 
most cases we assume that this phase relationship takes 
place, but do not introduce new denotations for V and V± 
for the sake of simplicity. In this case, (T,9 = — ^arg), in 
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correspondence with the formula l|42|) . the variances V± 
depend only on the difference between phases AO and we 
arrive at 

V = l + 2(n- | (aio 2 ) |), (46) 

fl = 2 Re((a 2 ) e^*^ ) - 2 | <a+a 2 ) | . (47) 

For the NOPO without additional polarization mixing 
(0+02) = (a 2 ) = (a 2 .) = and hence the case of the 
symmetric variances V+ = V- is realized. The phase- 
locked NOPO generally has non symmetric uncertainty 
region. However, the variances V+ and V_ become equal 
for the special case of 9 2 — 0i — $ x = § , when the insep- 
arability condition reads as V < 1, (V = V- = V + ). We 
note, that the relative phase <I> X plays an important role 
in specification of the entanglement. 

Entanglement in the case of unitary time-evolution 

So far, we have considered mainly steady-state regime 
of generation, including the effects of dissipation and 
cavity-induced feedback. However, in order to better 
understand the pecularities of the entanglement for the 
system under consideration, it would be interesting and 
desirable to study the case of a short interaction time 
t « 7 _1 , when the dissipation in the cavity is unessen- 
tial. In our analysis we shall assume that all processes 
are spontaneous, i.e. both subharmonic modes are ini- 
tially in the vacuum state. Then, the system's state in 
the absence of losses is generated from the vacuum state 
by the unitary transformation 

I* (t)) = u (t) 10), |0) 2 = e-***** |0)j |0) 2 , (48) 

which for x ~> gives two-mode squeezed vacuum state. 
Using formula (|43[) . where n = (u^ 1 (t)a^'aiu(t)) and 
{ai(t)a 2 (t)) = (u^ 1 (t)aia 2 u(t)) , and also considering the 
case of zero-detunings Ai = A 2 = A = 0, after a long 
but straightforward algebra we get the final result for the 
variance V = \(V+ + V-) for two different operational 
regimes: e < x an< 4 £ > X- The variance for the range e < 
X of a comparatively weak scaled pump field e = fc-E/73 
reads as 

e 2 

V(t) = 1 - -5 [cos(2/it) - 1] 
M 

- - sin(2jut) cos (£0) , (49) 
A 4 

where fi = \J x 2 ~ £ 2 - Thus, we observe a periodic evo- 
lution of V in this regime typical for the linear coupling. 
The level of squeezing of the two modes is periodically 
repeated. The behavior of the two-mode variance also 
significantly depends on the phase matching condition, 
so that we may tune the phase sum to maximize the 




FIG. 1: Unitary evolution of the variance V (t) versus the 
scaled dimensionless interaction time xt f° r the range of e < \ 
and provided that cos (SO) = 1. The parameters are: e/x = 
0.1 (curve 1), e/x = 0.4 (curve 2) and e/x = 0.7 (curve 3). 

entanglement. We further choose for illustrations the 
phase sum corresponding to cos (E0) = 1 for both oper- 
ational regimes. The dependence of V versus the scaled 
time-interval is shown in Fig. ^ where the three curves 
correspond to three different choices of the ratio e/x- 
Common to all curves is that the variance is nonclassi- 
cal and squeezed at least at the points of its minima £ m i n , 
which can be obtained by the formula ctg (2/it m j n ) = e//i. 
In all cases, the maximal degree of two-mode squeezing 
V m \ n = V (t) |t=t min — 0.5 is achieved in the limit e — ► x- 
For the curves on the Fig. ^ the time-intervals corre- 
sponding to the minimal values of V m i n equal to xi m in — 
0.74 + 1.0057rfc(curve 1), x^min ^ 0.63 + 1. 09 Ink (curve 
2), ximin ^ 0.57 + 1.47rft(curve 3), (k = 0, 1, 2, ...). 

If the opposite inequality holds, e > x 5 then the non- 
linear parametric interaction becomes dominant over the 
linear coupling and the variance is given by the following 
formula 

e 2 

V(t) = 1 ~~\~~ 77 [cosh (2ryi) - 1] 
r\ A 

- - sinh (2rft) cos (£0) , (50) 
77 

where r\ = y/ ' e 2 — X 2 ■ 

For x ~ * and cos (£0) = 1, from this formula we 
arrive at a well-known result, V (t) = exp (— 2et), for 
two-mode squeezed state. This shows that in the limit of 
infinite squeezing the corresponding state approaches to a 
simultaneous eigenstate of X\ — X 2 and Yi+Y 2 , and thus 
becomes equivalent to the EPR state. Fig. [21 shows the 
behavior of V when the system operates in the regime e > 
X- This figure clearly shows that as the interaction time 
increases, the variance decreases and reaches its mini- 
mum. Then the squeeze variance exponentially increases 
with the growth of the interaction time. For the data 
in Fig. |21 we obtain that the time intervals for which the 
variance reaches the minima are xi m i n ~ 0.48(curve 1), 




0,0 ^ , , , , , , , , , , 

0,0 0,2 0,4 0,6 0,8 vf 1,0 



FIG. 2: Unitary time evolution of the variance V (t) for 
the range of e > x aim provided that cos (SB) = 1. The 
parameters are: e/x = 1-1 (curve 1), e/x = 2 (curve 2) and 
e/x — 3 (curve 3). 

X*min — 0.38 (curve 2), x*min — 0.31(curve 3). It is easy 
to check that these points of minima can be found by the 
formula ctgh(2i]t mul ) = e/r), provided that cos = 1. 

We also conclude that the variance squeezes up to a 
certain interaction time, if e > X- It should be noted that 
although the time evolution of the variances are quite 
different for each of the operational regimes, the minimal 
values of the variance are described by the formula which 
is the same for both regimes 

Vmin = V(t)| t=tmta = -£-. (51) 

e + x 

With increasing e/x, in the regime e > x, the minimal 
value of the variance decreases as V min ~ x/e <C 1, which 
means that perfect squeezing takes place in the limit of 
infinite pump field. This result is not at all trivial for 
the system considered, even in the absence of dissipation 
and cavity-induced feedback, because the insertion of po- 
larization mixer usually destroys the two-mode squeezing 
produced by nondegenerate parametric down-conversion. 
The reason is that two-mode squeezed vacuum state is a 
supperposition of two-photon Fock states \n} 1 |n) 2 and 
the polarization mixer destroys the Fock states having 
the same number of photons, i.e., n\ — n-i — n. The de- 
tailed analysis of this problem can be found, for example, 
in [3 



1,0 n 

V 
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0,8- 
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0,6- 
0,5- 
0,4- 

1 2 3 4 5 

t/t th 

FIG. 3: Minimized variance V versus dimensionless amplitude 
of the pump field e/eth = kE/j3^/ (x — |A|) 2 + -y 2 for both 
operational regimes. The parameters are: xll — 0.1, A/7 = 
10 (curve 1), xll — 0.5, A/7 = 3 (curve 2) and xH ~ 
0.5, A/7 = 1 (curve 3). 



It is easy to check that in the limit x — > the vari- 
ance coincides with the analogous one for the ordinary 
NOPO, V = V- = V+ = 1 - e/(e + ^7 2 + A2 )- We 
see that the minimal variance remains less than unity for 
all values of pump intensity and is a monotonically de- 
creasing function of e 2 . For all parameters the maximal 
degree of two- mode squeezing V ~ 0, 5 is achieved within 
the threshold range. It is also easy to check that this 
expression is well-defined for all values of e 2 , including 
the vicinity of the threshold. One should keep in mind, 
however, that the linear approach used does not describe 
the threshold range where the level of quantum fluctua- 
tions increases. As a consequence, some matrix elements 
of (|29|l . (|39() . I|40|) increase infinitely in the vicinity of 
threshold. Nevertheless, it follows from l|52l) and further 
results l|54 |) -(|57 |) . that such infinite terms are cancelled in 
the variance V, as well as in VL , V+ for both operational 
regimes. This is not surprising, since such cancellation 
of infinities in the quadrature amplitude variances takes 
place for the ordinary NOPO also. 




Sub-threshold regime 

Using formulas (|29fl. ||I5J, after some algebra, we ob- 
tain the minimal variance in the following form 



V = 1 



e ( eS 2 - v/ 7 2 5 4 + A 2 (S 2 - 2 X 2 )' 
Si - 4A 2 x 2 



(52) 



One of the differences between the squeezing effects 
of the ordinary and self-phase locked NOPO is that the 
variances VL and V+ for the ordinary NOPO are equal 
to each other, while for the self-phase locked NOPO they 
are in general different. The values of the non symmetric 
variances are expressed through R according to formula 
(|42H . This quantity can be calculated with the help of 
the formula l|47|l . The result is found to be 
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R 



5-4 _ 4 A 2 % 2 



( 7 2 + (A- x ) 2 ) ( 7 2 + (x + A) 2 ) 



7 2 S 4 + A 2 {S 2 - 2 X 2 ) 



2e 



(53) 



As we see, the final expressions Q52|. (|53[1 below the 
threshold are rather unwieldy. We show the correspond- 
ing numerical results on the Figs. for illustration. 



Above-threshold regime 

Performing calculations for each of the cases A > 
and A < 0, we find the variance V in the above-threshold 
regime in the following form 



1,0 



v 3 - 1 

4 Vl + (s 2 - 



X 



4J/7 2 



4 I A 



(54) 



The results (|52f) and (|54|l for both operational regimes are 
summarized in Fig. [21 where the variance V is plotted 
as a function of the amplitude of the pump field. One 
can immediately grasp from the figure that the sum of 
variances V = ■k(V+ + VL) remains less than unity for all 
nonzero values of the pump field, provided that \ < A|. 
This shows the nonseparability of the generated state. 
The maximal degree of entanglement is achieved in the 
vicinity of threshold, V ~ 0, 5, if xl |A <C 1. In far above 
the threshold, E >> Eth, V increases with mean photon 
numbers of the modes and reaches the asymptotic value 
y = 3/4 + x/4|A|. It should also be mentioned that the 
result H54|) is expressed through the scaled pump field 
amplitude e = fc-E/73 and hence depends on coupling 
constants k and x- 

As the last step in connecting two-mode squeezing to 
observables of CV entanglement, we report the expres- 
sions of the non symmetric variances calculated with the 
help of the formulas |@2J|, Pol l and (gTJ. Upon evaluat- 
ing all required expectation values, we obtain 



R 



sgn{A) \A\-x 



1 



I A | V l + {£ 2_ £ih)h 2 



(55) 



and hence 



V + = 



3 | A | +x + sgn (A) cos (Afl) (| A | - x ) 

4|A| 
1 + sgn (A) cos (A0) 



V- = 



4 v /l + ( £ 2 -^)/7 2 ' 



3 | A | +x — sgn (A) cos (A0) (| A 

4i^n 

1 - sgn (A) cos (A0) 



(56) 



(57) 



v.,v 



0,9 
0,8 
0,7 
0,6 
0,5 
0,4 



(2) 



FIG. 4: The variances V+ (curve 1) and V- (curve 2) versus 
the scaled pump field amplitude e/eth for the parameters: 
Xh = 0-5, A/7 = 3. 



For the case of A0 = 0, when the variances are maximally 
different, the results are reduced to 



V- = 



1 



X 

2A : 



V+ = l- 



1 



2 v /l + (e 2 -^)/ 7 2 



(58) 



for A < 0. The case of A > is obtained from (|58|l by 
exchanging V+ ~ + V- and V- — * V+. These results for 
both operational regimes are depicted on Fig. 0] Some 
features are immediately evident. First of all one can 
see that both variances are minimal in the critical range 
but show quite different dependences on the ratio e/eth 
above the threshold. Attentive reader may ask about de- 
pendence of the results obtained on parametric coupling 
constant k. We note in this connection that the squeezed 
variances are expressed through the scaled pump field 
amplitude e = kE/j 3 and hence depend in general on 
both coupling constants k and x- 

In terms of demonstrating the CV strong EPR entan- 
glement, one has to apply another criterion V+V- < \. 
For the case of the symmetric uncertainties (V+ = V- = 
V), the product of the variances V 2 > \ and hence the 
strong EPR entanglement can not be realized. In the 
general case, the product of the variances reads as 



V+V- = V 2 -R 2 cos 2 



(A0) 



(59) 



It seems that V+V- lies below 1/4 at least for the relative 
phase A0 = ±7rrn, (m = 1,2, ...), and in the vicinity of 
threshold, where V ~ 0,5. However, for such selection of 
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v + v » 
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0,6- \ 
0,5- ^ 

0,4 -I . , . , . , . , . , 

1 2 3 4 , 5 

th 

FIG. 5: The product of the variances V+ and V- versus the 
scaled pump field amplitude e/e t h for the parameters xll = 
0.1, A/7 = 10 (curve 1) and xll = °- 5 ' A /7 = 1 (curve 2). 

the phases we arrive at 

v+v _ = iAi+x( 2 - 1 y (60) 

4 I A I V 

It is easy to check that the product of the variances ex- 
ceeds \ even in the vicinity of the threshold. This quan- 
tity for both operational regimes is illustrated in Fig. [5] 
It should be noted again that a detailed analysis of this 
problem, must include more accurate consideration of 
quantum fluctuations in the critical ranges. 

CONCLUSION 

Our work demonstrates the possibility of creation of 
CV entangled light-states with well-localized phases. We 
show that such, so-called, entangled self-phase locked 
states of light can be generated in NOPO recently re- 
alized in the experiment [|| . This device is based on the 
type-II phase-matched down-conversion and additional 
phase localizing mechanisms stipulated by the intracavity 
waveplate. The novelty is that this device provides high 
level of phase coherence between the subharmonics in 
contrast to what happens in the case of ordinary NOPO, 
where the phase diffussion takes place. This development 
paves the way towards the generation of bright CV entan- 
gled light beams with well-localized phases. It looks like 
that this scheme involving phase locking may be poten- 
tially useful for precise interferometric measurements and 
quantum communications, because it combines quantum 
entanglement and stability of type - II phase matching 
with effective suppression of phase noise. The price one 
has to pay for these advantages is the small aggrava- 
tion of the degree of CV entanglement in comparison 
with the case of ordinary NOPO. The quantum theory 
of self-phase locked NOPO has been developed in linear 
treatment of quantum fluctuations for both below- and 



above-threshold regimes of generation. We have studied 
the CV entanglement as two-mode squeezing and have 
shown that entanglement is present in the entire range 
of pump intensities. In all cases the maximal degree of 
two-mode squeezing V ~ 0, 5 is achieved in the vicin- 
ity of the threshold. It has also been shown that the 
amount of the entanglement can be controlled via the 
phase difference $ x . The other peculiarities of the sys- 
tem of interest have been established for the case of uni- 
tary dynamics. One of these concerns the presence of 
two operational regimes generating two-mode squeezing. 
If the linear coupling between subharmonics dominates 
over the parametric down-conversion, e < Xi we have 
observed a periodic evolution of the squeezed variance. 
Maximal degree of squeezing has been V ~ 0,5 in this 
regime. If the parametric interaction becomes dominant, 
e > X, the more high degree of two-mode squeezing can 
be obtained, Vmi n = < 0.5, but up to certain inter- 
action time. 

In our analysis we have not investigated all possible 
quantum effects of self-phase locked parametric dynam- 
ics. In particular, we have noted that the system con- 
sidered displays different types of quantum correlations, 
but we have not analyzed their connection with all pos- 
sible kinds of the entanglement. Consideration of quan- 
tum fluctuations in the near threshold operational range 
of self-phase locked NOPO also deserves special atten- 
tion for more accurate identification of CV entanglement. 
These topics are currently being explored and will be the 
subject of forthcoming work. 
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